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Abstract 

We present the first experimental study of the electric field distributions En of a three- 
dimensional (3D) microwave chaotic rough billiard with the translational symmetry. The transla- 
tional symmetry means that the cross-section of the billiard is invariant under translation along z 
direction. The 3D electric field distributions were measured up to the level number N = 489. In this 
way the experimental spatial correlation functions Cat iP (x, s) cx (EV^x + ^s)E^ (x — ^s)) were 
found and compared with the theoretical ones. The experimental results for higher two-dimensional 
level number iVj_ appeared to be in good agreement with the theoretical predictions. 
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In this paper we present the first experimental investigation of electric field distributions 
of the chaotic 3D microwave rough billiard with the translational symmetry. Due to experi- 
mental difficulties there are very few experimental studies devoted to 3D chaotic microwave 



cavities 



In a pioneering experiment Deus et al. |l( have been measured eigen- 
frequencies of the 3D chaotic (irregular) microwave cavity in order to confirm that their 
distribution displays behavior characteristic for classically chaotic quantum systems, viz., 
the Wigner distribution. Three-dimensional chaotic cavities as well as properties of random 
electromagnetic vector field have been also scarcely studied theoretically [y, |7|, [8|. 

In general, there is no analogy between quantum billiards and electromagnetic cavities in 
three dimensions. However, for 3D cavities with the translational symmetry the classification 
of the modes into transverse electric (TE) and transverse magnetic (TM) is possible. The 
TM modes are especially important because they allow for the simulation of 2D quantum 
billiards on cross-sectional planes of 3D cavities. Furthermore, we show in this paper that 
the distributions of the electric field of TM modes of the 3D chaotic rough cavity can be 
experimentally measured. 

In the experiment we used 3D cavity with the translational symmetry in the shape of 
a rough half-circle (Fig. [1]) with the height /?, = 60 mm. The cavity was made of polished 
aluminium. 

We suppose that the direction of the translational symmetry of the cavity is along the 
z-axis. The boundary conditions at z = and z = h demand that the z dependence of 
the ^-component of the electric and magnetic fields En, p (x) and -Btv jP (x) of TM modes be 
in the form E N>p (yL) = E NjP (x,y,z) = A NjP tp NtP (x,y)f p (z), where f p (z) = cos(pirz/h), p = 
0, 1, 2 . . ., An )P is the normalization constant and -Bjv iP (x) = 0. The functional dependence 
of £jv,p(x) on the plane cross section coordinates is denoted by the amplitude ipN tP (x,y) = 
Ejsr !P (x,y). The amplitude ipN, P {x,y) satisfies the Helmholtz equation 

(A ± + k 2 N>p )^ N>p (x,y) = 0, (1) 

where Aj_ is two-dimensional Laplacian operator and k^, P = {^ — (pn/h) 2 ) 1 / 2 is the effective 
wave vector. The wave vector = 27rz/jv/c, where z/jy is the resonance frequency of the 
level N and c is the speed of light in the vacuum. The equation (1) is equivalent to the 
Schrodinger equation (in units h = 1) describing a particle of mass m = 1/2 with the kinetic 
energy kj^ in an external potential V = (pn/h) 2 Therefore, microwave 3D cavities 
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FIG. 1: Upper panel: Sketch of the chaotic half-circular 3D microwave rough billiard in the xy 
plane. Dimensions are given in cm. The cavity sidewalls are marked by 1 and 2 (see text). Squared 
wave functions \ipN,p{Rc, 9)\ 2 were evaluated on a half-circle of fixed radius R c = 9.25 cm. Billiard's 
rough boundary T is marked with the bold line. Lower panel: White circle marks the position of 
the hole drilled in the upper wall of the cavity. This hole was used to introduce the perturber 
inside the cavity in order to measure the z-component of the electric field distributions £V iP (x). 



with the translational symmetry simulate on the cross-sectional planes quantum billiards 
with the external potential (pir/h) 2 . In this way microwave cavities can be effectively used 
beyond the standard 2D frequency limit (the case p = 0) jlOj in simulation of quantum 
systems. The amplitude ipN,p(%,y) fulfills Dirichlet boundary conditions on the sidewalls of 
the billiard. Therefore, throughout the text the amplitudes iPn{x, y) are also often called the 
wave functions ipN(x,y). It is important to note that the full electric field £^jv,p( x ) satisfies 
additionally Neumann boundary conditions at the top and the bottom of the cavity. 

Because of the relatively low quality factor of the cavit y (Q ~ 4000) the value of the level 



number iV was evaluated from the Balian-Bloch formula 



ill 
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where k is the wave vector, V = (9.43 ± 0.01) -10 4 m 3 is the volume of the cavity and 
J s ^ = 0.932 m ±0.005 m is the surface curvature averaged over the surface of the cavity. 
The measurements allowed us for the first time to evaluate the spatial correlation function 



12] 
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C NtP (x, s) = (Ejv iP (x + -s)£^p(x - -s)), (3) 



where the local average (■ • •) is defined as follows 

(l^, P (x)| 2 ) = \ N ] \E N ^ + s)\ 2 d n s. (4) 

lA n J-A/2 

The 3D cavity sidewalls are made of 2 segments (see Fig. [1]). The rough seg- 
ment 1 is described on the cross-sectional planes by the radius function R(9) = Ro + 
J2m=2 a m sm(m9 + (f) m ), where the mean radius i?o=10.0 cm, M = 20, a m and m are uni- 
formly distributed on [0.084,0.091] cm and [0,27r], respectively, and < 9 < tt. (Here, for 
the convenience, the polar coordinates r and 9 are used instead of the Cartesian ones x and 

y) 

The surface roughness of a billiard on the cross-sectional planes is characterized by the 

function k(9) = {dR/ d9) / Rq. For our billiard we have the angle average k = ((k 2 (9)) e ) 1 ^ 2 ~ 

0.400. In such a billiard the classical dynamics is diffusive in orbital momentum due to 

collisions with the rough boundary because k is much above the chaos border k c = M -5 / 2 = 

0.00056 [14]. The roughness parameter k determines also other properties of the billiard [ijj] 

on the cross-sectional planes. The amplitudes ipN,p{ r , 9) are localized for the two-dimensional 

level number N ± < N e = 1/128& 4 , where N ± = £k 2 Np -£k NjP . A = (1.572±0.002)-10~ 2 m 2 

and P = 0.537 m ±0.001 m are the cross-sectional plane area and its perimeter, respectively. 

Because of a large value of the roughness parameter k the localization border lies very low, 

iV e ~ 1. The border of Breit-Wigner regime is Nw = M 2 /A8k 2 ~ 52. It means that between 

I L 

N e < N± < N w Wigner ergodicity [15j ought to be observed and for N± > N w Shnirelman 
ergodicity should emerge. 

To measure the amplitudes ^jy^fc 9) of the 3D electric field distributions we used a 



very effective method described in 16f| . It is based on the perturbation technique 171 

nnn 

and preparation of the "trial functions" [18J, H9|, 12DJ • In the perturbation method a small 
perturber is introduced inside the cavity to alter its resonant frequencies and in this way to 
evaluate the squared wave functions \ipN, P (R c , 9)\ 2 (see Fig. [p. The perturber (4.0 mm in 
length and 0.3 mm in diameter, oriented in z-direction) was moved by the stepper motor 



FIG. 2: The reconstructed wave function ^460, o( r >$) of the chaotic half-circular microwave rough 
billiard. The amplitudes have been converted into a grey scale with white corresponding to large 
positive and black corresponding to large negative values, respectively. Dimensions of the billiard 
are given in cm. In the figure the z dependence of the electric field distribution E^Q^{r, 9, z)\ x= o oc 
ip460,o(r,0)\ x=0 fo(z) is also shown. 

via the Kevlar line hidden in the groove (0.4 mm wide, 1.0 mm deep) made in the cavity's 
bottom wall along the half-circle R c . The measurements were performed at 0.36 mm steps 
along a half-circle with fixed radius R c = 9.25 cm. 

In order to find the dependence of the electric field distributions E^^fe) on the z coor- 
dinate and to estimate the wave vector fc 3 = pn/h we measured the electric field inside the 
3D cavity along the z-axis. Also in this case the perturber (4.5 mm in length and 0.3 mm in 
diameter) was attached to the Kevlar line and moved by the stepper motor. The perturber 
entered and exited the cavity by small holes (0.4 mm) drilled in the upper and the bottom 
walls of the cavity. The both holes were located at the position: r = 9.11 cm, 9 = 0.47 
radians. 

Using the method of the "trial wave function" we were able to reconstruct 75 experi- 
mental wave functions ipN, P (j',9), which belonged to TM modes of the rough half-circular 
3D billiard with the level number N between 2 and 489. The range of corresponding eigen- 
frequencies was from u 2 ~ 2.47 GHz to z/ 48 9 ~ 11.99 GHz. The remaining wave functions 
belonging to TM modes, from the range N = 2 — 489, were not reconstructed because of 
near-degeneration of the neighboring eigenfrequencies or due to the problems with the mea- 
surements of \ipN,p(Rc, 9)\ 2 along a half-circle coinciding for its significant part with one of 
the nodal lines of ipN, P (r : 9). 
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FIG. 3: The reconstructed wave function ^463,4 ( r > °f the chaotic half-circular microwave rough 
billiard. The z dependence of the electric field distribution #463,4 (r, 0, z)\ x= o oc "0463,4 (r, 0)\ x= ofn(z) 
is also shown. 



In Fig. [2] and Fig. [3] we show two examples of reconstructed wave functions ^460,0 i r ,9) 
and ^463,4 (r, respectively. The character of the wave functions predominantly depends 
on the effective wave vector k^,p- It is seen that the wave function ip463,4 (r, 8) in Fig. [3] is 
more regular than the one presented in Fig. [2] in spite of having the larger level number N. 

In order to check ergodicity of the billiard's wave functions ipN, P (R c ,9), especially close 
to the ergodicity borders, one should use some additional measures such as e.g., calculation 
of the structures of their energy surfaces 14j. For this reason we extracted wave function 
amplitudes C^' p ^ = (n,l\N,p) in the basis n, I of a half-circular billiard with radius r max , 
where n = 1, 2, 3 . . . enumerates the zeros of the Bessel functions and I = 1, 2, 3 . . . is the 
angular quantum number. As expected, close to the border of the regimes of Breit-Wigner 
and Shnirelman ergodicity the wave function ipmw{r, 9) (N± = 65) was found to be extended 
homogeneously over the whole energy surface [13] (figure not shown here). In contrary, the 
wave function ip^^r, 0), N± = 16, which lies closer to the localization boarder, displays 
the tendency to localization in n, I basis (figure not shown here). 

The measurement of 3D electric field distributions En iP (x) allowed us for the first time 
to find the experimental spatial correlation function Cjv iP (x, s). It is easy to show [21] that 
for the 3D chaotic cavity with the translational symmetry the spatial correlation function 
should have the following form: 

Cjv, P (x, |s|) = Cjv, P (|s|) = J (k NjP s xy ) cos(p7rs z /h), (5) 

where |s| = (s 2 xy + si) 1 / 2 . For the cross-sectional planes z = const the correlation function 
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FIG. 4: Panels (a)-(c) show the experimental correlation function C46o,o( x > l s l) calculated at x = 
(-2.75 cm, 4.35 cm, cm) for the three projection angles = 0, vr/4, and tt/2, respectively. 
Experimental correlation function 6460,0 ( x > l s l) (full line) is compared with the theoretical one 
(dashed line). 



Cjv,p(l s l) ~ Jo(kN,pS X y) is reduced to the well known result of Berry [21] for chaotic 2D wave 
functions described by a random superposition of plane waves. 

In Fig. HJa)-(c) we show a representative example of the experimental correlation function 
C460,o(x, |s|) (N± = 65) calculated at x = (-2.75 cm, 4.35 cm, cm) for the three different 
projection angles = 0, 7r/4, and tt/2, respectively, where = arcsin(s z /|s|). The local 
average (• • ■) required for the evaluation of Cjv, p (x, |s|) (see the formulas (3) and (4)) was 
calculated on the cross-sectional plane xy in the range A/2 = 2n/kN, P - The experimental 
correlation functions C*46o,o( x ) l s l) ar e compared in Fig. H]with the theoretical ones. In all 
cases we find good agreement with the theoretical predictions given by the formula (5). 



Small discrepancies observed in Fig. IU(a) for 
the system and were theoretically studied in 



s| > 1 can be connected with the finiteness of 



22] 



7 




'0 2 4 

|s| [cm] 



FIG. 5: Panels (a)-(c) show the experimental correlation function C463,4(x, |s|) calculated at x = 
(-2.75 cm, 4.35 cm, cm) for the three projection angles = 0, vr/4, and tt/2, respectively. 
Experimental correlation function 6*463,4 (x, |s|) (full line) is compared with the theoretical one 
(dashed line). 

Fig. Eta)-(c) shows the experimental correlation function C463,4(x, |s|) (N± = 16) cal- 
culated at x = (-2.75 cm, 4.35 cm, cm) for the three different projection angles = 0, 
7r/4, and 7r/2, respectively. The experimental correlation functions C463,4(|s|) are compared 
in Fig. [5] with the theoretical ones. In Fig. [5] (a), even for small |s|, we find a significant 
departure of the experimental correlation function from the theoretical prediction, which 
clearly suggests that the wave function ipiQ3^(r, 6) is not chaotic. Also in Fig. [5](b)-(c) the 
experimental correlation functions C463,4(|s|) show for larger |s| significant deviations from 
the theoretical ones. The discrepancies between the correlation function 6463,4 (x, |s| = z) 
for the z-component of the electric field distribution and the theoretical prediction in Fig. 
El^c) arise mainly due to the procedure of averaging of the correlation function Cat, p (x, |s|) 
in z-direction, which was taken over the period of the cosine function. 



S 



In summary, we measured the wave functions of the chaotic 3D rough microwave billiard 
with the translational symmetry up to the level number N = 489. For the first time the ex- 
perimental correlation function CV jP (x, s) was estimated and compared with the theoretical 
prediction. For the states with higher N± we find, especially for small values of the parame- 
ter |s|, good agreement with the theoretical predictions, which show that the wave functions 
are chaotic. For the states with lower N± significant discrepancies between experimental 
and theoretical results are observed. 
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